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Abstract 

The experimental finding of an ubiquitous kink in the nodal direction of angle-resolved photoe- 
mission spectroscopies of superconducting cuprates has been reproduced theoretically. Our model 
is built upon the Migdal-Eliashberg theory for the electron self-energy within the phonon-coupling 
scenario. Following this perturbative approach, a numerical evaluation of the bare band disper- 
sion energy in terms of the electron-phonon coupling parameter A allows a unified description of 
the nodal-kink effect. Our study reveals that distinction between A and the technically defined 
mass-enhancement parameter A* is relevant for the quantitative description of data, as well as for 
a meaningful interpretation of previous studies. 

A remarkable agreement between theory and experiment has been achieved for different samples 
and at different doping levels. The full energy spectrum is covered in the case of LSCO, Bi2212 
and overdoped Y123. In the case of underdoped Y123, the model applies to the low energy region 
(close to Fermi level). 

PACS numbers: 74.25. Gz, 74.25.Kc, 74.72.-h, 79.60.-i 
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I. INTRODUCTION 



High resolution angle-resolved photoemission spectroscopy (ARPES) is nowadays con- 
sidered as one of the most powerful methods for obtaining detailed information about the 
electronic structure of atoms, molecules, solids and surfaces.- With the advent of improved 
resolution, both in energy and momentum, ARPES has provided key information on the elec- 
tronic structure of high temperature superconductors, including the band structure, Fermi 
surface, superconducting gap, and pseudogap.-^ In this context, its large impact on the de- 
velopment of many-body theories stems from the fact that this technique provides a means 
of evaluating the so-called electron self-energy, S. A wide and comprehensive document 
on the more relevant aspects of the photoemission spectroscopy is found in Ref. Q. As a 
description of the spectroscopic techniques based on the detection of photoemitted electrons 
is beyond the scope of this paper, below we will only summarize some aspects that will be 
relevant in the course of our discussion. 

As regards the photoemission process, although it formally measures a complicated non- 
linear response function, it is helpful to notice that the analysis of the optical excitation of 
the electron in the bulk greatly simplifies within the "sudden-approximation" In brief, 
this means that the photoemission process is supposed to occur suddenly, with no post- 
collisional interaction between the photoelectron and the system left behind.- In particular, 
it is assumed that the excited state of the sample (created by the ejection of the photo- 
electron) does not relax in the time it takes for the photo-electron to reach the detector.- 
It can be shown that within the sudden approximation using Fermi's Golden Rule for the 
transition rate, the measured photocurrent density is basically proportional to the spectral 
function of the occupied electronic states in the solid, i.e.: Jk oc A^(E). Eventually, and 
validated by whether or not the spectra can be understood in terms of well defined peaks 
representing poles in the spectral function, one may connect A^{E) to quasiparticle Green's 
function G(k, E) = 1/(E — Ttk(E) — with Ttk(E) defining the self-energy and Ek the bare 
band dispersion. In fact, A(k,u) = — (l/7r)ImCr(k, u + i0 + ). Beyond the sudden approxi- 
mation, one would have to take into account the screening of the photoelectron by the rest 
of system, and the photoemission process could be described by the generalized golden rule 
formula, i.e, a three-particle correlation function.- For our purposes, it is important mention 
the evidence that the sudden approximation is justified for the cuprate superconductors even 
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at low photon energies.- 1 ^ In the end, the suitability of the approximations invoked, will be 
justified by the agreement between the theory and the experimental observations for the 
wide set of data. 

Interactions involving a low-energy excitation appear as a sudden change in the electron 
energy dispersion near the Fermi level (Ep), known as kink.— This feature has been observed 
in various CSC and at different doping levels both along the nodal^^^^^^ 1 ^ 1 ^ 1 ^ 1 ^ 
and antinodal^^J^Ji'JJi directions, revealing diverse and controversial behaviors commonly 
interpreted in terms of the coupling of electrons to a certain kind of excitation. Recall that 
the antinodal direction denotes the (n, 0) region in the Brillouin zone where the d-wave 
superconducting gap has a maximum. This fact extraordinarily complicates the theoretical 
interpretation of the kink, due to the anisotropic character beyond the s-channel.- More 
advantageous is the nodal direction that corresponds to the (0,0)-(7r,7r) direction in the 
Brillouin zone, where the c?-wave superconducting gap is zero. Multiple measurements have 
been realized along this direction in several CSC both in the normal and superconducting 
(SC) state. Such measurements show a kink in a similar energy scale (in the range of 48-78 
meV) and are present over an entire doping range, and for temperatures well below and 
above T c . 

To the moment, there is no consensus on the origin and behavior of the kink in the CSC 
and its influence on the bosonic coupling mechanism that leads to SC state. Even more, 
there is no model allowing to reproduce the kink effect in different materials and/or different 
doping levels. For this reason, a theoretical description of the kink effect in various CSC and 
at different doping levels is quite desirable. In this article, with the aim of obtaining a correct 
dressed electron band dispersion relation which can reproduce a wide number of experimental 
spectroscopies, we suggest a model based on the Migdal-Eliashberg (ME) approach for the 
numerical determination of the electron self-energy S. As a main result, we emphasize the 
importance of considering the proper distinction between the mass-enhancement parameter 
A* and the electron-phonon coupling parameter A. 

The essence of the Eliashberg theory is a perturbative scheme that allows to deal with 
strong electron-phonon coupling effects. It relies on the Migdal approach for the theory of 
metals, that is posed in terms of thermal Green's functions G for systems of many interacting 
particles described by a hamiltonian H within the Fermi liquid picture.— 1 ^ Eventually, 
one is lead to the definition G~ l = Gq 1 — E, where Gq stands for the non-interacting 
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electrons, whereas E is the above referred self-energy function which is a measure of the 
perturbation introduced to the bare Hamiltonian by the interactions. Remarkably, the 
averaging procedures for obtaining E happen to be expressed in terms of experimentally 
accessible data (u 2 F) convoluted with well known mathematical functions (R), i.e.: E oc 
a 2 F * R. a 2 F will be derived from inelastic neutron scattering experiments and R defined 
in terms of the so-called digamma functions. In this work, the comparative study of a wide 
set of experimental results has allowed an empirical extension of the commonly considered 
physical scenario when analyzing ARPES data. 

The paper is organized as follows. First, in SecHIl we will give details about the relevant 
quantities and theoretical treatment (ME approach) to be used. The physical interpretation 
of the underlying approximations is also focused. In Sec lIII[ the analysis of the ARPES 
data according to our proposal is done. Comparison with previous published material will 
be emphasized. Finally, Sec JIVI is devoted to discuss our results. The relevance of the 
electron-phonon coupling mechanism for the interpretation of the nodal-kink of ARPES 
experiments in CSC will be concluded. 



II. THE ELECTRON SELF-ENERGY E 

In ARPES the dressed electronic dispersion relation is denoted by E^.— This quantity 
characterizes the charge carriers as quasiparticles that are formed when the electrons are 
dressed with excitations. E^ has been commonly related to the bare band dispersion 
through the real part of the electron self-energy by E^ = + ReH(Ek). Within the Eliash- 
berg theory, the electron-phonon interaction (EPI) self-energy may be obtained from the 
real part of the expression^ 



V \2 2ttT J v \2 2vrT J J ' 

(1) 

valid for the whole range of temperatures T, frequencies v, and energies u. Here, ip(z) are 
the so-called digamma functions with complex argument and a 2 F{v) defines the important 
EPI spectral density which measures the effectiveness of the phonons of frequency v in 
the scattering of electrons from any state to any other state on the Fermi surface. More 
specifically, being interested in the nodal direction ARPES experiments, which are not 
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influenced by the anisotropy of the superconducting gap, we will refer to a (non-directional) 
isotropic quasiparticle spectral density, defined as the double average over the Fermi surface 
of the spectral density a 2 F(k, k', v>), i.e.: 

o?F{v) = ^ E I sL> I' - <-Me k )5(e kr ) (2) 

kk'J 

where, g kk , = [h/2MuL i ] 1 ' 2 (k|e k , k ■ W|k') defines the EPI-matrix element for electron 
scattering from k to k' with a phonon of frequency f k _ k > (j is a branch index). M stands for 
the ion mass, V is the crystal potential, e k , k is the polarization vector, and N(0) = J2k $( £ k) 
represents the single-spin electronic density of states at the Fermi surface. As usual S(x) 
denotes the Dirac's delta function evaluated at x' = 0. Owing to the intrinsic complexity for 
evaluating a 2 F(u) for a given material from first principles, in this work we have adopted a 
more empirical point of view, which relies on auxiliary experimental data. 

Taking into consideration the inherent existence of phonons in the CSC, and for simplicity, 
we have chosen the model by Islam & Island for the extraction of the EPI spectral function 
a 2 F(u) through the phonon density states obtained from inelastic neutron scattering in 
LSCO,— Bi2212,— and Y123.— The corresponding spectral densities are shown in Fig. [TJ in 
comparison with the spectral densities obtained by Shiina and Nakamura,— and Gonnelli et 
aL— We emphasize that very similar results are found under the use of any of these densities. 

At this point, we should comment that other choices of the spectral density are possible, 
but have been left aside. For instance, the absence of the magnetic- resonance mode in 
LSCO and its appearance only below T c in some CSC are not consistent with the idea that 
the nodal kink has a magnetic origin.- Along the same line, the reported absence of the 
magnetic-resonance mode in Bi2212 at a doping level of 0.23 without the disappearance of 
the superconducting stated (T c m 55K) reveals that the magnetic-resonance mode cannot 
be directly related to the spectral density involved in Eq. Therefore, we have preferred 
the spectral densities oriented to phonons. 

An important ingredient of our model is the commonly used dimensionless electron- 
phonon coupling parameter, defined in terms of the EPI spectral function by A = 
2 J °° dv a 2 F [y) j v . This quantity is not to be straightforwardly identified with the mass- 
enhancement parameter A*. As it will be shown below, coincidence is only warranted under 
certain limits. A further relevant feature to mention is that Ifi 1 ^! 2 and as a consequence 
A are inversely proportional to the number of charge carriers contributed by each atom of 
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FIG. 1: (color online) EPI spectral density a 2 F(v) for LSCO, Y123 and Bi2212. The black solid 



lines correspond to the method of Ref 
the methods of Ref. 
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The other lines are shown for comparison employing 



25l in all materials (dashed lines), and from Ref. 



2J in Bi2212 (dotted lines). 



the crystal to the bosonic coupling mode. Therefore, an increase in the doping level, which 
causes an increment in the hole concentration of the CuO<i plane must be reflected in a 
reduction of A as we will see in the analysis of the kink structure. Furthermore, recalling the 
outstanding feature of the theory of metals, that Ifi^/) 2 vanishes linearly with |k — k'| when 
|k — k'| <C k F ,— one would expect a linear disappearance of the coupling effect that gives 
rise to the nodal kink in the vicinity of the Fermi surface. On the other hand, inspired by 
recent results on the universality of the nodal Fermi velocity vp < (at low energies) in certain 
cuprates, a prominent role of this quantity is also expected. 



A. Green's function formalism 



Some words are due, considering Eq.([T]) and the relation between the dressed and bare 
energies. Within the arguments customarily used for analyzing ARPES data, one finds 
the low energy approximation (close to the Fermi surface) Eh ~ £fc/(l + A), with A the 
electron-phonon coupling parameter defined above. However, one should recall that such 
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expression is an asymptotic form of the more correct Ek ~ £fc/(l + A*) with A* the true 
mass-enhancement parameter and A an appropriate limiting value. As a central result of 
our work, distinction between them is essential for the overall description of available data. 
Let us show how this arises. 

Following the thermal Green's function formalism, we assume that the electron-phonon 
interaction is introduced by the relation 

G' 1 {k, iuj n ) = G^ 1 (k, iu n ) - E(k, iuj n ) (3) 

with Gq 1 related to the bare electron energy and iu n standing for the "imaginary Matsubara 
frequencies" .— 

We recall that, technically, the bare electron band energy Ek is determined by the poles 
of the Green function G(k,iu n ), or the zeros of G~ 1 (k,iu n ) at the iu n frequencies.— On 
the other hand, it is known that additional dynamical information is contained in the an- 
alytic continuation G(k,u + i0 + ) to points just above the real frequency axis, known as 
the "retarded" Green's function. One is therefore led to continue the electronic self-energy 
E(/c, lo + i0 + ) analytically by T,(kuj + i0 + ) = to) + iT: 2 (k, lj). Now, suppose that a pole 
occurs near to = so, one gets 

G- 1 (k,uj + m + ) =w-e fc -Ei(A;,u;) -iE 2 (k,u) (4) 
~ u (l - ^£^[ =o ) - [en + S^fc, 0)] - iE 3 (fc, w) . 

Then, the pole of G occurs at a frequency cuq given by uq = Ek — i/2rk, with r^ 1 = 
-2(l-^S 1 )- 1 S 2 (A;, J E fc ) and 

^(l-^EO-^efc + EiCM)] • (5) 

Here A^ = — ^w^iL=o * s the technically defined mass-enhancement parameter.— 

We want to emphasize that replacement of A* by A in Eq. ([S]) is only warranted for 
st cites k qXi or very close to the Fermi surface when the ME approach for the self-energy 
[Eq.flTJ] has been employed at the low temperature limit (T — > 0). Again, owing to the 
difficulties for evaluating Ei beyond the Migdal approximation, our position in this paper 
has been to obtain A* through the systematic evaluation of ARPES data, as shown in 
Sec JIIBI As a further detail, related to the limitations introduced by the use of the electron- 
phonon coupling parameter, one should mention that, long before the advent of the high 



T c superconductivity, Ashcroft and Wilkins showed that in some simple metals the single 
parameter A is insufficient to determine a number of thermodynamic properties such as the 
specific heat.— 

In the following, we will introduce the simplest correction possible, for dealing with the 
correct mass-enhancement parameter. Noteworthily, it will be shown that a phenomenolog- 
ical linear relation, i.e.: A* ~ 5X suffices for the interpolation of experiments and numerical 
data. The physical interpretation of the parameter 5 will be discussed within the following 
section flHBj) . Just from the technical side, we want to comment that the mathematical 
material within this section has been developed using the convention of positive energy 
bands relative to the Fermi surface. It is apparent that the contrary selection can also 
be done, changing signs within intermediate expressions but the same final results, i.e.: 
Ek — (1 + cLEi)^ 1 [ek — 0)], with Sk defining negative energy bands. 

B. Phenomenological dispersion relation 

Let us make some final remarks before moving onto the application of the above ideas 
to the analysis of the ARPES data. Let us start by recalling that the bare electron band 
energy Ek is not directly available from the experiments. Instead, the electron momentum 
dispersion curve E%(k — kp) may be measured. However, it has been noted that the relation 
between the dressed and bare energies is a central property as related to the kink structure. 
In fact, based on the commonly used equation E). = ej t + Re{T J {Ei c )} one can consider that 
implicitly depends on Ek through the boson coupling parameter A.— ^ This fact, along with 
some ansatz for the ARPES "bare" dispersion allows to obtain A as a unique constrained 
parameter that better fits the observed kink structure. However, the indiscriminate proposal 
of dispersion relations could considerably under or overestimate the average renormalization 
because implicit approximations are used as indicated in the previous section. Furthermore, 
we stress that, the use of the same ansatz on different cuprates or even for a definite material 
with slight variations in the doping level is not warranted. In this work, with the aim of 
finding a widespread renormalization function that describes the behavior of the kink in 
different CSC and for different doping levels, we have carried out an exhaustive study on 
the incidence of the EPI coupling parameter in the appearance of kink-dispersion. Based on 
an interpolation scheme between the numerical behavior of the relation Ek = £fc+i?e{£(-Efc)} 
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and the experimental data, as a central result, we have encountered that the practical totality 
of data are accurately reproduced by a universal dispersion relation of the kind 

k-k F = ^(l- — x] , (6) 

v F< V ^log / 

with v F< the Fermi velocity at low-energies, c^i og the so-called logarithmic frequency as 
introduced by Allen,— and w the (only) free parameter required for incorporating the specific 
renormalization for each superconductor. Note that the constant frequency co>i og (properly 
defining the corresponding spectral densities a 2 F(v) for each A) has been introduced only for 
scaling purposes, just with the aim of reducing the scattering of numerical values in dealing 
with different samples. In this sense, w is a mere mathematical instrument. Thus, our 
numerical program is as follows: (i) v F< is determined from the momentum dispersion curve 
in the ARPES measurements [~ 1.4 eVA - 2.2 eVA], (ii) u\ og is evaluated from the spectral 
densities shown in FigJT]through the definition u;i og = exp{(2/\) J °°ln(z/)[a 2 F(z/)/i/](ii/} (we 
get u{f g co ~ 16.1455 meV, c^ 23 ~ 35.5900 meV and wg 2212 ~ 33.8984 meV respectively), 
(iii) Ek is numerically determined from the relation E k = Ek+Re{'E(Ek)}, and (iv) correlation 
is established between theory and experiment by the application of Eq.((6]). 

From the physical point of view, our empirical ansatz [Eq.flEJ)] may interpreted as follows. 
Let us assume that the involved quantities are not far from their values at the Fermi level, 
and start with e k replaced by E k — Tii(E k ), i.e. 

E k - H x (E k ) = (k — k F ) v F< (1 - 5XY 1 , (7) 

where the definition 5 = vj/oj\ og has been used. Now, let us take derivatives respect to E k 
and evaluate for E k — > 0. One gets 

d(k - k F ) 



1 + A* 



v F< (l-5\y\ (8) 

E k =0 



dE k 

and recalling that v F< is obtained as the slope of the lower part of the momentum dispersion 
curve,— this equation leads to 1 + A* = (1 — dX)^ 1 . Thus, a physical interpretation of the fit 
parameter 8 is obtained, i.e.: 5 = (A*/ A) / (1 + A*). To the lowest order, the dimensionless pa- 
rameter S is basically the ratio between the defined mass-enhancement and phonon-coupling 
parameters 5 ~ A*/A. Outstandingly, it will be shown that this fact reassembles the differ- 
ences obtained by tight-binding Hamiltonian models 3 ^ 3 ^ and the "density-functional" band 
theories^ 2 - 13 ^. Recall that, in principle, the density-functional theory 3 ^ gives a correct ground- 
state energy, but the bands do not necessarily fit the quasi-particle band structure used to 
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FIG. 2: (color online) The renormalized energy as a function of momentum k — kp for 
La2- x Sr x CuO/± with various doping levels between 0.03 (right curve) and about 0.30 (left curve), 
all measured at a temperature of 20.ff along the nodal direction. Experimental data are taken from 
Ref. |lCj. The theoretical curves are labeled according to the best fit values for A. 



describe low-lying excitations. As it will be seen below, predictions from both types of 
models may be reconciled appealing to the differences between A and A*. 



III. ANALYSIS OF ARPES DATA 



Below, we present the application of our theoretical analysis for a wide set of experimental 
curves available in the literature. The main facts are shown in figures [2]- HJ arid summarized 
in tables imi 

a. Results for LSCO.- In Fig. [2] we show the results found in La 2 - x Sr x CuOi covering 
the doping range (0 < x < 0.3). Remarkably, within this range, the physical properties span 
over the insulating, superconducting, and overdoped non-superconducting metal behavior. 
Superconducting transition temperatures T c in the interval of 30-40K have been observed 
by Bednorz and Miiller— and others^ 1 ^. For the application of Eq.Q, here, we have 
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TABLE I: Values of the boson-coupling parameter A and the corresponding mass-enhancement 
parameter A* obtained from the analysis of ARPES data at several doping levels of La2- x Sr x CuO^. 



A* has been 


obtained by means of Eq. ([ 


3) to the lowest order A* ? 


a 5 A (in this case 5 = 


0.185). The 


predicted superconducting transition temperatures T c are also shown. Our results are 


presented in 


contrast with other models available in 
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1.90 


0.35 


37.83 




0.18 


1.80 


0.33 


36.19 




0.22 


1.50 


0.28 


30.47 




0.30 


1.30 


0.24 






0.1-0.2 [34J 6 


2-2.5 




30 - 40 




[25] 6 


1.78 




40.6 




0.15 [32j c 


1 - 1.32 


0.14 - 0.22 






0.22 [321 c 


0.75 - 0.99 


0.14 - 0.20 





a Wc allow a margin of error in A of ~ ±0.3 as related to the numerical interpolation procedure between 



theory and experiment. 



b The A 
c In Ref. 



/allies reported in that reference were obtained so as to fit T c at the indicated values. 
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the electronic structure of LSCO has been calculated employing a generalized gradient approx- 



imation to density functional theory and used to determine A. 

considered vp < = 2eV-A as related to the experimental results of Refs. On the 

other hand, the best fit of the whole set of experimental data has been obtained for 5 = 0.185, 
and the derived A* values are shown in the Table [B For comparison, recall that values of 
"A = 2 — 2.5" in the range 0.2 > x > 0.1 were reported in the Ref. |3J by Weber. In that 



case they were obtained within the framework of the nonorthogonal-tight-binding theory 
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of lattice dynamics, based on the energy band results of Mattheiss^- and corresponding to 
the range of the EPI coupling parameter A predicted in our model. It must be emphasized 



that in the case of Ref. 
LSCO. Moderate discrepancies 
the analysis of Refs. 



34 and 



3J, A was obtained in agreement with the observed T c values in 
Detween the predictions of our phenomenological model and 
may be ascribed to some uncertainty in the experimental 
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spectral densities. As regards the critical temperatures, we have calculated them based on 
the celebrated McMillan's equation,^ T c = (w Iofl /1.2)escp[-1.04(l + A)/(A - fj*(l + 0.62A))] 
(see tabled]). In all the calculations, the Coulomb's pseudopotencial was given a typical 
value, /i* = 0.13. It is essential to be aware that, there is no small parameter which 
enables a satisfactory perturbation theory to be constructed for the Coulomb interaction 
between electrons. Thus Coulomb contributions to the electron self-energy S cannot be 
reliably calculated.— Fortunately this is not a serious problem in superconductivity because 
a reasonable assumption is to consider that the large normal-state Coulomb effects contained 
in the Coulomb self-energy are already included in the bare band structure The remaining 
off-diagonal terms of the superconducting components of the Coulomb self-energy turn out 
to have only a small effect on superconductivity, which is treated phenomenologically.— 
One can see that the consideration of the electron-phonon interaction in LSCO strongly 
suggests that the high T c values can be caused by conventional electron-phonon coupling, in 
agreement with the conclusion of Weber.— 

From a different perspective, in a recent publication, Giustino and co-workers^ have 
calculated the electronic structure of LSCO employing a generalized gradient approximation 
to density functional theory (DFT). These authors have extracted the "A" parameter by 
measuring the gradients of both the theoretical and the experimental self-energy data within 
the low energy limit (0 — 50meV). Their procedure yields \ exp t = 1-00 — 1.32 for the 
optimally doped sample (x = 0.15) at 20K and \ exp t = 0.75 — 0.99 for the overdoped 
sample (as = 0.22) while the theoretical results A^ = 0.14 — 0.22 at optimal doping and 
Xth = 0.14 — 0.20 in the overdoped regime. It is concluded that theoretical values noticeably 
underestimate the experiments and, thus, that the electron-phonon interaction is unlikely 
to be relevant in LSCO. From our view, considering that the theoretical value is calculated 
from the gradient of the DFT self-energy, X t h is basically to be identified with A*, while X eX pt 
matches the electron-phonon coupling A involved in the standard Migdal formalism analysis 
of experiments. 
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FIG. 3: (color online) Same as in Fig. [2] but in samples of Bi2Sr2CaCu20s+ x with (left to 
right)::E = 0.12 (underdoped), x = 0.16 (optimally doped), and x = 0.21 (overdoped). The red 



17. The continuous lines 



diamonds and blue circles correspond to the experimental data of Ref. 
correspond to our theoretical curves. 

b. Results for BSCCO.- In Fig. [3] we display the results found in samples of Bi2212. 
Each case has been studied with temperatures both in the normal and superconducting 
state. The experimental data were taken from the work by Johnson et al— In this case^we 
have used vf < = 1.6eV-A as a value consistent with the experimental results of Refs. 
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171 The best fit with experimental data has been found for 5 = 0.354. Similarly to the 
case of LSCO, our analysis fits well the "A" values predicted by others and very different 
models^ 1 ^ 1 ^ (see table HT])- From our results, it is clear that the phonon coupling mode as a 
unique source for the behavior of the critical temperature in BSCCO is only reasonable for 
the underdoped case (x=0.12). Nevertheless, the full energy spectrum in the nodal direction 
and the consequent emergence of the kink-effect are reproduced. 

c. Results for YBCO.- Finally, in Fig. [4] we show the results found for Y123 samples. 
The experimental data were taken from recent work by Borisenko et alM To our knowledge, 
no more experimental evidence of kinks in the nodal direction is available for Y123. The 
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TABLE II: Same as tableUbut for the case of Bi2Sr2CaCu20$ +x .Hei:e, several values of A are not 
related to ARPES data, but just introduced for comparison. In this case, 5 = 0.354 



SC x 


Ref. 


A 


A* 


TJX) \K] 

C\ V L J 


Bi2212 0.12 


This" 


2.15 


0.76 


64.81 


0.16 




1.33 


0.47 


42.45 


0.21 




0.85 


0.30 


19.93 




[25f 


3.28 




85 




This 


3.28 


1.16 


81.66 




[26] 6 


3.34 


1.05 


93 




This 


3.34 


1.18 


82.40 


0.16 


[12j c 


~ 1.28 


~ 0.43 





a Wc allow a margin of error in A of ~ ±0.3 as related to the numerical interpolation procedure between 
theory and experiment. 

h Thc A "patties reported in that reference were obtained so as to fit T c at the indicated values. 

c In Ref. 1 12 two channels are defined for A. Ai = 0.43 ±0.02 corresponds to the "primary" channel (close to 

Fermi level) and is free from normalization effects. A2 is obtained from the Kramers-Kronig transformation 

and the experiment^ In this sense, and within the notation employed in the current work, we obtain, 

A* = Ai, and A ~ 0.85 + 0.43. 

value vf < = 1.63eV-A has been used for consistency with the experimental results reported 
by those authors. The best fit between the experimental data and our model has been found 
with the value 5 = 0.365. It must be noted that the appearance of a second kink in the 
underdoped case may not be allocated in the current model, and requires further theoretical 
considerations. Contrary to the properties of the overdoped case, that is captured by the 
theory in the full energy spectrum, the underdoped case is only reproduced at low energies 
(1st kink close to the Fermi level). One possibility for upgrading the theory would be 
to consider the existence of additional contributions to the electron self-energy £ which 
could be related to high orders of the phonon perturbation, vertex corrections, or inclusive 
Coulombian effects. On the other hand, one specific feature of Y123 is that, contrary 
to the other superconductor cuprates, the reservoir layers in these materials contain CuO 
chain layers which could be contributing significantly to the band energies of the in-plane 



14 



■100 



> 

CD 
LU 




■150 



-200 



-250 



0.12 0.06 



Exp. 
X = 0.65 



T=30 K 
x = 0.6 



Exp. 
^=0.50 



T=30 K 
x = 0.85 




0.0 0.12 0.06 0.0 0.12 
k-k p (A 1 ) 



0.06 0.0 



FIG. 4: (color online) Same as in Fig. [2] but in samples of Y Bd2Cuj > 0%j rX with (left to right): 
x = 0.4 (underdoped), x = 0.6 (underdopedV, and x = 0.85 (overdoped). The solid squares 



la . The lines correspond to our theoretical curves. 



correspond to the experimental data of Ref. 
All curves have been obtained at 30K. 

electronic structure. This hypotheses has been used with success by Cucolo et al.— for the 
interpretation of tunneling spectra, specific heat and the ultrasonic attenuation coefficient 
in both phases of Y123. More detailed studies on the electronic photoemission spectra of 
Y-based copper oxides are required. Unfortunately on the basis of the existing data, then, it 
is not possible to favor one or other of these possibilities. Moreover, in ARPES one should 
also care about the residual 3-dimensionality and its effect on photoemission data. Again, 
our model fits well the "A" values predicted by other works^ 1 ^ 1 ^ (see table UTTj) . However, 
is clear that the phonon coupling mode as the source for the critical temperature in YBCO 



samples is not reasonable. The same conclusion has already been obtained in the Refs. 
and 36. 



35 
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TABLE III: Same as tables U and [TT] but for Y Ba^Cu^O^x samples. In this case, 5 = 0.365 



SC x 


Ref. 


A 


A* 


Tc(X) [K] 


YBCO 0.4 


This a 


0.80 


0.29 


17.51 


0.6 




0.65 


0.24 


9.82 


0.85 




0.50 


0.18 


3.39 


- 


[25j 6 


~3.4 


- 


91 


— 


This 


3.45 


1.26 


84.19 




[35J 6 


~0.5 




~3 




[35J 6 


~1.3 




~30 




This 


1.30 


0.47 


36.43 




[36J C 




0.18 - 0.22 






This 


0.49 - 0.60 


0.18 - 0.22 


~3.0 - 6.6 



a We allow a margin of error in A of ~ ±0.3 as related to the numerical interpolation procedure between 
theory and experiment. 

h The A vphles reported in that reference were obtained so as to fit T c at the indicated values. 

36l the parameter "A" has been obtained from the spectral density a 2 F(k, v) employing the 



c In Ref. 



local density approximation to density functional theory (see text). This value correspond at the mass- 
enhancement parameter A* 



IV. CONCLUDING REMARKS 



In summary, we have introduced a model that allows to reproduce the appearance of the 
ubiquitous nodal kink for a wide set of ARPES experiments in cuprate superconductors. 
Our proposal is grounded on the Migdal-Eliashberg approach for the self-energy of quasi- 
particles within the electron-phonon coupling scenario. The main issue is the proposal of 
a linear dispersion relation for the bare band energy, i.e.: = (k — kp)vp < (l — <5A) -1 . 6, 
the only free parameter of the theory is a universal property for each family of cuprates. It 
has been interpreted as the relation between the mass-enhancement A* and electron-phonon 
coupling A parameters. 

An excellent agreement between the theory and the available collection of experiments 
is achieved. Our results support the idea that the phonon coupling mechanism is the main 



1(3 



cause of the kink effect and the so-called universal nodal Fermi velocity, although its effect 
in the appearance of the superconducting state and the high critical temperatures is not 
clear yet. 

For decades, a well-known controversy has arisen on the role of the parameter "A" whose 
values noticeably scatter among different model calculations. As a central result, our pro- 
posal re-ensembles the "A" values obtained from different models and, as a first approx- 
imation, solves the controversy through the relation A* = 5X. We emphasize that the 
phenomenological parameter 5 (obtained through the analysis of a wide collection of data) 
has allowed to go beyond the conventional Migdal-Eliashberg approach analysis of restricted 
sets of experiments. 



36 



Our model is directly supported by the "A" values obtained in Refs. 
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25 



26 



32 



34 



35, 



When inserted into the celebrated McMillan's formula these values indicate that the 
electron-phonon coupling behind is not necessarily the only mechanism responsible for the 
superconducting properties of all the cuprates. In fact, the critical temperatures that one 
can calculate are only at reasonable levels for the case of LSCO (tabled!) and for underdoped 
BSCCO (Table [II]) where the phonon mechanism dominates. 
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